Conductivity engineering of graphene by defect formation 
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Transport measurements have revealed several exotic electronic properties of graphene. The 
possibility to influence the electronic structure and hence control the conductivity by adsorption or 
doping with adatoms is crucial in view of electronics applications. Here, we show that in contrast to 
expectation, the conductivity of graphene increases with increasing concentration of vacancy defects, 
by more than one order of magnitude. We obtain a pronounced enhancement of the conductivity 
after insertion of defects by both quantum mechanical transport calculations as well as experimental 
studies of carbon nano-sheets. Our finding is attributed to the defect induced mid-gap states, 
which create a region exhibiting metallic behavior around the vacancy defects. The modification of 
the conductivity of graphene by the implementation of stable defects is crucial for the creation of 
electronic junctions in graphene-based electronics devices. 



PACS numbers: 73.61.Wp, 72.80.-r 

I. INTRODUCTION 

The electronic properties of graphene have been of con- 
siderable interest lately, with several repor ts of unique 
structural and electronic properties! 1 * 2 * 3 * 4 ^! the next 
phase of research on graphene, it is most likely that we 
will witness an increased focus on the tailoring of elec- 
tronic properties by means of defects, adatoms and geo- 
metrical confinement. The possibility to open band gaps, 
induce gap states by defects and the ability to control the 
conductivity by chemical and physical means is hence 
a crucial step in order to establish graphene as a com- 
petitive material for future electronics. It has in fact 
already been shown that a bandgap can be opened in 
graphene, and that it can be tuned as a function of car- 
rier concentration. 5 The high sensitivity of graphene due 
to changes in carrier concentration is related to its low 
Johnson noise. This property was e.g. used in the cre- 
ation of a gas detector capable to detect the adsorption 
of single atomP Though adsorbed gases sensibly change 
the bandstructure of graphene, they tend to be unstable 



on the graphene surface. Molecular doping of NO2 and 
N2O2 induce p-type behaviour in graphene P The maxi- 
mum resistance levels are lower for higher concentration 
of doping. Another road to engineer the bandstructure of 
graphene consists in the addition of more stable defects 
using wet chemistry methods. This is attractive in view 
of electronics applications since such defects are physi- 
cally more stable than gases attached to the graphene 
layer. In fact, it has been shown recently that chemically 
introduced vacancy defects change the electronic struc- 
ture of grapheneP 

From the XAS analysis on a graphene layer which was 
chemically treated in the same manner as in this work, 
we have observed a spectroscopic signature that could 
be attributed to defect sites in the graphene planes (i.e. 
of a* character) P Though the distribution of defects in 
graphene is certainly complex and a matter of current 
studies, these states of a character related to a defect 
site could only be oriented in the graphene plane, and 
it has been analyzed as being the result of the presence 
of a vacancy. Hence we start our study by considering 
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how vacancy defects influence the electronic structure of 
graphene. For this reason we have considered both single 
vacancy and di- vacancy in the ir bonded graphene struc- 
ture through a tight-binding model with nearest neigh- 
bour interactions. In viewing the graphene as consisting 
of two sub-lattices, sub-lattice A and B, it is known that 
the effect of an impurity in the A sub-lattice is mani- 
fested in the B sub-lattice!^ In particular, a vacancy in 
the 7r-band, i.e. absence of tt orbitals at the impurity 
site generates mid-gap states for the atoms in the B sub- 
lattice located in a neighborhood around the vacancyP 
These mid-gap states arise due to symmetry breaking 
that removes the equivalent Dirac points in the two sub- 
lattices. 

Defects always decrease the mobility of sample com- 
pared to its defect-free version, while they may increase 
the number of carriers. The latter effect is especially 
important here and leads to a strong conductivity rise, 
since for pure graphene there are no carriers at fermi en- 
ergy. We explore the relevance of defects for controlling 
the conducting properties, using both theoretical and ex- 
perimental approach. In our theoretical study we inves- 
tigate the effects on the local density of states (LDOS) 
from single vacancies and di-vacancies by means of real 
space Green functions (GFs). In order to understand the 
influence of the electronic structure on the conductivity, 
we calculate the evolution of the LDOS as a function of 
distance from the defect position. From there we pro- 
ceed with ab-initio density functional calculations of the 
electronic structure and conductivity. The experimental 
studies consist of point contact conductivity measure- 
ments on free-standing carbon nanosheets (CNS). The 
CNS consists of few-layer graphene with mono- or double 
layers at the edges (Fig. [lj. The samples in this study 
were grown on a conducting substrate that serves as back 
contact, greatly facilitating the experimental setup. De- 
fects were introduced by wet chemical treatment, in ac- 
cordance with previous workP 



II. MODELLING OF DEFECTS AND 
PREDICTION OF PROPERTIES 

A. Modified local density of states 

The local density of states (LDOS) of a single graphene 
sheet, both acid treated and untreated, is modeled in 
real space by the dressed graphene Green function (GF) 
G r (r, r'; u) = Gg(r - r'; u) + G r (r - r ; W )T(u,)G&(r - 
r; uj); here we discuss the results for an impurity located 
at the position i"o = {r^r^}. The dressed GF is thus 
given in terms of the homogeneous GF Gg(r — r';ui) 
with the diagonal and off-diagonal elements given around 
the Fermi level Ep = by g\(v — r';ui) w — Jo(kp\r — 
r'|)[2wlog(£>/M) + *7rTw]/(47i>w|) and g r AB {v - r';w) « 
— Jo(kp\r— r'|)(2£)+i7rw)/(47Tj(w|,), respectively, whereas 
the T-matrix is given by 




FIG. 1: Scanning electron micrograph of CNS a) on the 
reference sample at 52° tilt and b) on the acid treated sample 
at 52° tilt, both grown on tungsten substrates c) Schematic 
image of in-situ sharpened probe needle in contact with a CNS 
d)SEM image of the contact between the omniprobe needle(l) 
and a free standing carbon nano sheet (2). 



Here, U = (Ua,Ub) is the impurity potential matrix, 
where Ua (Ub) is the scattering potentials at (rg), 
whereas Jq{x) is a Bessel function of the first kind, vp 
is the Fermi speed, kp is the Fermi wave vector, and 
p is the graphene planar density and is related to the 
energy cut-off D ~ 5 — 10 eV through Airpvp = D 2 . 
Although this approximation does not describe the mod- 
ified LDOS quantitatively correct, it gives a sufficient 
qualitative account. First we consider a single impurity 
located at ya- As we are considering a vacancy, we let 
Ua — * oo. The perturbation to the LDOS in the A sub- 
lattice then becomes SNa(tc,lj) — —^sSG r A (r,r;uj)/TT 
uj Jq (kp\r — r a\) / (An pvp) , showing that the Dirac points 
remain in the A sub-lattice. The perturbed LDOS in the 
B sub-lattice, on the other hand, becomes 
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T(w) = [1 - UG^O.uj)]- 1 ^. 



(1) 



which is clearly diverging at the Fermi level. This di- 
vergence signifies the existence of mid-gap states for the 
carbon atoms in the B sub-lattice. Using more sophisti- 
cated approximations, which removes the unphysical sin- 
gularity, show that the mid-gap states are indeed present 
in the B sub-lattice PD 



B. Calculations of electronic structure and 
resistivity 

We performed ab-initio density functional calculations 
of the graphene electronic structure in presence of sin- 
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gle vacancies and di-vacancies. The calculations em- 
ployed the projector augmented wave method, and the 
plane wave cut-off energy was set to 750 eV. The gener- 
alized gradient approximation (GGA) was used for the 
treatment of exchange-correlation functional and a lat- 
eral 8x8x1 supcrccll of graphcnc with a divacancy was 
considered for the calculations. This yields an effective 
vacancy concentration of 0.0156. The geometry of the 
supercell in presence of the defect was optimized using 
Hellmann-Feynman forces with a tolerance of 0.01 eV/ 
A. A T-centered 3x3x1 set of k-points was used in the 
Monkhorst-Pack scheme with a Gaussian broadening of 
0.2 eV. 

In the calculations of the graphene conductivity we 
used an all-electron scalar-relativistic version of the tight- 
binding linear muffin-tin orbital (TB-LMTO) method™ 1 
within the local spin-density approximation 11 (LSD A) to 
the density functional theory.^ The valence basis con- 
sisted of s-, p-, and d-orbitals. The atomic sphere ap- 
proximation (ASA) was employed here, which required 
an additional empty sphere to be present in the graphcnc 
layer per two carbon atoms and other empty spheres be- 
tween graphene layers with a spatial distribution similar 
to the one already used successfully for graphite P^l The 
integrations over the 2D BZ were performed on a uni- 
form mesh of about 5000 ky points. Such high density of 
ky-mesh is crucial for the accurate description of states 
near the Dirac point. 

The random distribution of defects in the system is 
described by means of the coherent potential approxi- 
mation. (CPA^SEIl This single-site approximation has 
been shown to be successful in predicting properties of 
systems with random disorder, thus averaged over all 
possible combinations of impurity distributions. The 
CPA can deal with systems with arbitrary low concen- 
tration of impurities while its numerical requirements re- 
main modest compared to the supercell method for dis- 
ordered systems. The CPA has been used in other stud- 
ies of graphene, where it was argued that CPA gives very 
good result s for t he description of the physical properties 
of graphen c! 16 l 17 [ The question concerning applicability 
of the CPA near Dirac point was studied by Skrypnyk, 
Lokte\ff^. They have noted that only renormalized ap- 
proaches as CPA can be efficient close to the van Hove 
singularities in the spectrum. Furthermore, it was found 
that in the vicinity of the Dirac point a correction to the 
CPA self-energy is needed. We study DOS and resistance 
only for systems with a finite concentration of impurities, 
which naturally induces a shift of the Fermi level away 
from the Dirac point. Hence, CPA is applicable for the 
present study.™ The conductance was calculated from 
the Kubo linear response theory. 



III. EXPERIMENTAL VERIFICATION 



A. Graphene model system 

The results from theory are predictions for an ide- 
alized material which is not in contact with any sur- 
roundings. In order to experimentally verify the predic- 
tions, we have used well-characterized carbon nanoshccts 
( Figs, [lja) and [ljb)). These CNS were grown by 
plasma-assisted chemical vapor deposition (CVD) on a 
conducting substrate as a model system for graphene 
The CNS, consist of 1/im graphene sheets protruding 
from a 20 nm layer of graphite. The estimated thickness 
is 2-3 graphene layers at the top and about 7 monolayers 
in the vicinity of the substrate. The thus supported CNS 
is allowing both wet chemical treatment^ and point con- 
tact measurements without the adverse effects associated 
with removal of the CNS from the growth substrate and 
the deposition and contacting required for the conductiv- 
ity study (schematics is shown in Fig. [ljc))!^ Compared 
to graphene adhered and contacted i.e. on top of Si02, 
this set-up is likely to reduce the influence of substrate 
charges and impurities from graphene processing which 
might shade the effect of defects on the conductivity. 



B. Sample preparation 

The defects were introduced in the CNS by a treat- 
ment with concentrated aqueous HC1 as described,^ the 
CNS were exposed to 35% concentrated aqueous HC1 at 
90° C for 3 hours. The HC1 was removed from the flasks 
using a pipette and the samples were treated with deion- 
ized water for 10 h at 90°C, followed by a thorough rinse 
with RT dcionized water and drying in air at 150°C for 3 
min. In order to assess the effect of this acid treatment, a 
reference sample was prepared using deionized water in- 
stead of the aqueous HC1 then followed by the same wash 
- rinse - dry protocol. Water, in contrast to aqueous HC1, 
cannot induce permanent defects. The samples are in the 
following referred to as " acid-treated" and " reference" , 
respectively. All samples, includning the as-deposited, 
were handled in air. 



C. DC measurement setup 

The CNS samples were electrically contacted with 
an Omniprobe nano-manipulator inside a FEI Strata 
DB235 focused ion beam/scanning electron microscope 
(FIB/SEM). The Current- Voltage (I-V) measurements of 
the CNS were performed by using a Keithley 6430 source 
meter. The tungsten tip of the manipulator is in-situ ion 
polished in the FIB providing a better definition of the 
tip/graphene contact area as well as a reduction of the 
oxygen on the W surface. 

In our nano-manipulation experiment, the CNS are 
contacted with a sharp tip at room temperature. The 
probe needle is manoeuvred in the x, y and z-directions 
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until it contacts a CNS. All movements are carefully mon- 
itored in real time by continuously scanning the elec- 
tron beam in the FIB/SEM. The tip/CNS contact in 
Fig. [ljd) is established in the upper part of the CNS (for 
schematic representation of tip and carbon nano sheets, 
see Fig. |ljc)) . The tip/CNS contact can be precisely 
observed through a slight bending of the CNS in the 
realtime scanning electron micrograph (SEM) images. 
The voltage was scanned from -1 to IV in 30-40 seconds 
with no applied gate voltage. A reference SEM image 
was recorded before and after each measurement. When 
the tip/CNS contact is not sufficiently strong, occasional 
drift and tip vibrations led to a (I-V) temporary detach- 
ment of the tip seen as a current drop in the curves. 
Such I-V curves were not considered in the further evalu- 
ation. In this measurement geometry CNS are available 
in abundance with a surface density of 10 s cm" 2 . There- 
fore, the conductivity of a large quantity (> 400) of CNS 
could be assessed and consequently, the data can be ex- 
tracted with a high statistical precision. Three to five 
measurements are taken at the same location. The mag- 
netic field inside the FIB/SEM chamber is found to be 
approximately 66mT in Ultra High Resolution Mode. 
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IV. RESULT AND DISSCUSSION 
A. Influence of Vacancy defect on Density of States 

In our first approach to understand the electronic 
structure in functionalized graphene, the LDOS of a sin- 
gle graphene sheet, both defect free and defect contain- 
ing, is modeled by means of real space GFs. Consider a 
single impurity located at ta in sub-lattice A. Although 
the LDOS is slightly modified in sub-lattice A, the Dirac 
point remains. On the other hand, in sub-lattice B, the 
effect from the vacancy is more pronou nce d in that mid- 
gap states arise from the carbon atomsP^The emerging 
mid-gap state suggests that the B sub-lattice becomes 
metallic in presence of the vacancy in the A sub-lattice. 
Generalizing the above consideration to the case of a 
di-vacancy, with two adjacent vacancies of which one is 
in each sub-lattice, shows the emergence of the mid-gap 
states in both sub-lattices for atoms in a neighborhood 
around the di-vacancy. This conclusion indicates that 
graphene becomes metallic in a neighborhood of the di- 
vacancy. The length scale associated with the mctallic- 
ity around the di-vacancy is of the order of a few lattice 
parameters, Fig. [2] which is due to the spatial power 
law decay of the influence from the impurity. In order 
to substantiate the model considerations above we have 
performed density functional calculations using the Vi- 
enna Ab-initio Simulation Package for a graphene sheet 
with a di-vacancy! 22 1 23 1 The value of the DOS at the Fermi 
level (Ep), which signifies the existence of a metallic con- 
ducting component, is shown in the Fig. [3] A relatively 
large peak arising from the p z orbital develops at the 
Fermi level for an edge atom around the di-vacancy. If 



FIG. 2: (Color online) Real space distribution of the cor- 
rection SN(r,E F ) = -%J2i=AB SGi(r;E F )/Tv to the LDOS 
(arb. units) at the Fermi level in presence of the di-vacancy 
located at va = and rs = o(l, v3)/2 (large red dots). The 
coordinates are given in units of the lattice parameter a. 



one moves away from the defect site, the DOS at Ep 
decreases, but the metallic component extends over sev- 
eral lattice sites around the defect (more than 15 A away 
from the defect site). In the inset of Fig. 3, we show the 
DOS at the Fermi level for both sublattices, suggesting 
an oscillating behaviour. A closer inspection (main body 
of Fig. 3) reveals that this is simply due to that the two 
different sublattices in the graphene sheet show different 
decay rates away from the defect center. Our ab-initio 
calculations are in agreement with the model calculations 
discussed above, that is, vacancies introduce a metallic 
component in the electronic structure, and this metal- 
licity, extending to several lattice constants, should be 
sufficient to enhance the conductivity of the graphene. 



B. Influence of Vacancy defect on resistivity 

The central aspect of this paper is the conductance 
of graphene as a function of vacancy defects, which we 
describe in the following ab-initio calculations. For sim- 
plicity we only considered single C vacancies in these cal- 
culations, but in light of the discussion above, we expect 
the results to hold also for other vacancy geometries. 

The results of our calculations are shown in Fig. [4] It 
is obvious that as the concentration of vacancies grows, 
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Distance from hole edge (A) 

FIG. 3: Calculated Density of states (DOS) at the Fermi 
level for C atoms, as a function of distance from a di-vacancy 
of graphene. The values are shown for both sublattices. In 
the inset we show the same data, for both sublattices in one 
single graph. 



the OK temperature resistivity decreases with at least one 
order of magnitude compared to the material containing 
no vacancies. In the concentration regime of x < 0.005 
it is seen from Fig. [3] that the conductivity depends ex- 
ponentially with respect to defect concentration. This 
is actually a normal result for regular semiconductors^, 
and can be explained from a Fermi-Dirac distribution of 
electron states. For a larger concentration of vacancies 
(x > 0.005) it is found that the resistivity slightly in- 
creases (Fig. [2]), and appears to saturate close to a con- 
stant level for a defect concentration of 3-5 %. These 
results verify our conclusions drawn from the charac- 
ter of the electronic structure in presence of vacancies. 
The results in Fig. [4] should therefore be seen as a tran- 
sition from a semi-metallic regime (or rather a regime 
where the effects of Zitterbewegung dominate) with lim- 
ited conductivity^^ to a regime of highly conducting 
graphene. This transition is driven by defects via the 
creation of mid-gap states, which produce an extended 
regime of metallic character. Once this regime has been 
fully established, the addition of more defects produce 
scattering centers which reduce the conductivity, result- 
ing in a more conventional behavior. The connection be- 
tween defects and resistivity in graphene is very intricate. 
Similar changes in the density of states were observed in 
graphene nanoribbons, where states in the midgap ap- 
pear at interfaces P3 

Defects cause not only impurity states, but also accept 
or donate charge, thus causing a shift of the Fermi energy 
in the studied system. We have compared the Fermi lev- 
els for different impurity concentrations in order to quan- 
tify the influence of the Fermi level shift on conductivity. 
In order to distinguish between the contributions from 
the Fermi level shift and impurity scattering on the con- 
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FIG. 4: Calculated resistance of a single layer of graphene as a 
function of defect concentration. The red diamond represents 
a calculated resistance from a model consideration with a low 
amount of defects but a relatively large shift of the Fermi level 
(see text). 



ductivity, we have considered a model system as follows. 
The Green function of the model system is identical to 
that of a graphene with a very low concentration of de- 
fects (0.1 %), but the Fermi level corresponding to a much 
larger concentration of defects (metallic regime, i.e. 1 % 
defect) was used in the spirit of the rigid band approxi- 
mation. The system will further be denoted as graphene 
II. Therefore, the difference between the graphene II sys- 
tem and the original graphene consists only in a change of 
the Fermi energy. The calculated resistance of graphene 
II is shown as a (red) diamond in Fig. [4] Clearly the con- 
ductivity of graphene II is greatly enhanced compared to 
the original graphene with an impurity concentration of 
0.1%, although their only difference is 0.016 Ry change 
of the Fermi energy. 

In a metallic regime the residual resistance originates 
from impurity scattering and it is thus expected to be 
proportional to the impurity concentration. Graphene 
with x = 1% and the graphene II with x = 0.1% (Fig. [4} 
should differ mainly by the amount of impurity scatter- 
ing as compared to the different conductivity regime of 
graphene with x — 0.1%. When the ratio of the resis- 
tances between the two systems with the same Fermi en- 
ergy and different impurity concentration i.e. graphene 
at x = 1% and the graphene II at x — 0.1% is evaluated, 
it indeed approximates the expected proportionality for 
metals. Therefore the intial decrease of the resistivity 
curve in Fig. 4 is due to an increased metalic component 
of the electronic structure, with negligible influence from 
impurity scattering. At higher defect concentrations (0.5 
% and above) the effect of the impurity scattering be- 
comes important, and the balance between these two ef- 
fects causes the resistivity to grow with large defect con- 
centrations. 

In the experimental situation, the dangling bonds 
around the defects in graphene are most likely to be satu- 
rated with hydrogen, oxygen or hydroxide groups. In or- 
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FIG. 5: I-V characteristics in forward and reverse bias of the 
water treated (reference) sample. 



FIG. 6: Typical Current- Voltage characteristics in forward 
and reverse bias of a typical acid treated CNS. 



der to account for this, we have carried out density func- 
tional theory calculations saturating the dangling bonds 
around di- vacancies with H and/or O atoms. The LDOS 
of the carbon atoms around the defect at the Fermi level 
reaches, as in the calculations of vacancies with dangling 
bonds, finite values suggesting that metallicity persists 
upon addition of H and/or O to the carbon atoms in the 
neighborhood of the vacancies. 



C. IV characterization of Carbon nanosheets 

The presence of oxygen and absence of chlorine was 
verified by x-ray photoclcctron spectroscopy in model 
CNS system. In our XAS measurements, we observe in- 
planc, er* -type orbitals appearing after the acid treat- 
ment which indicates the presence of vacancies. In fact, 
when graphene is damaged for example by electr on irr adi- 
ation, both, single and di- vacancies were found pZHH B u t 
as shown above, the appearance of midgap states is valid 
for both single and di-vacancics and therefore, at this 
stage, our transport experiment is only sensitive to the 
presence of defects, but not to their detailed character. 
Fig. [5] shows three I-V measurements performed on the 
reference sample and Fig.[6]shows three I-V curves of acid 
treated (defect induced) sample, demonstrates that the 
I-V curves on individual CNS are reproducible. The I-V 
curves of Fig. [5] and Fig. [6] are symmetric and non- linear, 
where the reference sample shows currents in the range 
of 0.1-10 /xA and the acid treated sample in the range of 
0.01-1 mA at 1 V. The currents in the as-deposited sam- 
ple are about 2 times lower than in the reference sample. 
In order to evaluate the curves quantitatively, we have 
measured the resistance at bias voltages from -0.6V to 
0.6V. Every recorded I-V curve results in a number of re- 
sistance values taken at different bias voltages Vb ■ Both 
negative and positive voltages are considered for the eval- 
uation since no systematic asymmetry was seen in the I-V 
curves for opposite bias. The resistances are plotted in 
histograms for each voltage value (Fig. [7]). The scatter 
of the resistance data between the different CNS on the 



same sample are due to changes in size of the CNS and 
of the contact point of the needle. 

The measured data collected in resistance histograms 
were fitted with a log-normal distribution, which is used 
to determine the average resistance as well as the stan- 
dard error of the resistance at each voltage. These re- 
sistance values and their associated errors are shown in 
Fig. [8] as a function of bias voltage. The average resis- 
tance of the reference sample is about 50 times higher 
than the resistance of the acid treated sample. This 
change in resistance is much larger than the standard 
error of our measurement (Fig. [8]), and we can safely 
conclude that our described method of chemical treat- 
ment can be used to tune the conductance of CNS with 
more than one order of magnitude. Furthermore, from 
the experimental curves (Fig. [8| we observe for both ref- 
erence and acid treated samples a decrease of the resis- 
tance with increasing voltage. This could be either an 
intrinsic effect of the electronic structure of the CNS or 
could be related to local structural modifications of the 
nano sheets appearing at higher bias voltages. In fact, 
when we increase the voltage to about 4 V, we observe an 
increase in current by one order of magnitude followed by 
a drop in current close to the leakage current level of the 
setup which is about 50 pA. A subsequent SEM observa- 
tion demonstrates that in this case, several 10 nm large 
holes were created in the CNS. When the measurement 
was carried out between -1 V to 1 V, no visual changes 
could be noted on the CNS structures in the SEM im- 
ages and, as shown above, subsequent measurements on 
the same CNS show a reproducible resistivity. Therefore, 
the non-linearity of the I-V curves should be an intrinsic 
property of the CNS. 

From the SEM images, we measured the average sur- 
face area of the CNS in this study (Figs.Qa) and [jjb)). 
We obtain the resistivity to be 150 /i fim in the refer- 
ence sample and 3.5 /i f2m in the acid treated sample at 
0.3V. From a comparison with resistivity values reported 
in the literature, the resistivity of the reference sample 
corresponds well to the value found for nanocrystalline 
graphite with a crystallite size of 2-5nmP In fact, the 
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FIG. 7: Log histograms of resistances measured in typical 
CNS on the reference sample at a) 0.1V b) 0.5V and on the 
acid treated sample at c) 0.1V d) 0.5V. 



CNS studied in this work consist of structurally coherent 
domains or grains with a diameter between 2-5nm in very 
good correspondence to the above valued This similarity 
points to the CNS behaving as graphite concerning their 
resistivity. Deviations from a crystalline arrangement at 
C atoms in the grain boundaries in the as-grown CNS 
will lead to an increase in the resistivity as observed here. 
When a graphene sheet is exfoliated from a graphite crys- 
tal, the resistivity of this sheet when positioned onto a Si 
substrate will be 5 times lower, in the non-gated state, 
than in graphiteP Some consequences of the acid treat- 
ment on the CNS structure can be estimated from our 
observations. Structural changes with a size of several 
10 nm appear in the CNS after the acid treatment as ob- 
served in high resolution SEM images. Our conductivity 
measurements show that the acid treatment has modified 
the CNS structure. Such structural modification could 
be nucleated at step edges present in high density (the 
CNS thickness changes between 2-7 monolayers) and at 
the grain boundaries. One may speculate that a part of 
the about 50 times decrease of the resistivity upon acid 
treatment could be attributed to a thinning mechanism 
of the CNS, leaving on average more sheets exhibiting a 
graphene-like behavior. 

With this in mind, the decrease of the average resistiv- 
ity from 150 /i f2m to 3.5 /i f2m measured here, is by one 
order of magnitude higher than what would be expected 
from just a thinning of the CNS layers and moving from 
a graphite character to a graphene like conductivity. It 
should also be noted that the absolute value of the non- 
gated average resistivity is about 3.5 /i fim, which is as 
low as in measurements on low defect graphene P Still, 
the acid treated CNS will keep a similar grain size and 
grain boundary defects as the reference sample since the 
acid treatment is not likely to cure such defects. There- 
fore, we conclude that the huge change in resistivity upon 
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FIG. 8: Resistance (log scale) vs. applied voltage in approx. 
100 measurements (both forward and reverse bias are in- 
cluded) of reference (upper curve, solid line) and acid treated 
sample (lower curve, dotted line). The error bar is the stan- 
dard error from the log-normal distribution. 



acid treatment is largely driven by the creation of defects, 
and this experimental result is consistent with the the- 
oretical results reported here. The presence of vacancy 
defects is further stressed by our x-ray spectroscopic mea- 
surements which also show a very good agreement with 
the DOS calculated for vacancy defects in grapheneP 



V. CONCLUSION 

In conclusion we have demonstrated an increase in 
the conductivity of carbon nano-sheets and graphene by 
more than one order of magnitude by means of engi- 
neered introduction of defects. This behavior is drasti- 
cally different from the conductivity of defectuous carbon 
nano-tubesP3 The method to create the defects demon- 
strated in this paper is very simple and widely available, 
hence opening up a huge possibility of similar investi- 
gations on graphene in various geometries, on different 
substrates and in combination with other defects like ad- 
sorbed atomic or molecular species. 
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